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Abstract 

Following a previous work on the quantization of a massless scalar field in a space- 
time representing the head on collision of two plane waves which fucus into a Killing- 
Cauchy horizon, we compute the renormalized expectation value of the stress-energy 
tensor of the quantum field near that horizon in the physical state which corresponds to 
the Minkowski vacuum before the collision of the waves. It is found that for minimally 
coupled and conformally coupled scalar fields the respective stress-energy tensors are 
unbounded in the horizon. The specific form of the divergences suggests that when the 
semiclassical Einstein equations describing the backreaction of the quantum fields on 
the spacetime geometry are taken into account, the horizon will acquire a curvature 
singularity. Thus the Killing- Cauchy horizon which is known to be unstable under 
"generic" classical perturbations is also unstable by vacuum polarization. The calcu- 
lation is done following the point splitting regularization technique. The dynamical 
colliding wave spacetime has four quite distinct spacetime regions, namely, one flat re- 
gion, two single plane wave regions, and one interaction region. Exact mode solutions 
of the quantum field equation cannot be found exactly, but the blueshift suffered by 
the initial modes in the plane wave and interaction regions makes the use of the WKB 
expansion a suitable method of solution. To ensure the correct regularization of the 
stress-energy tensor, the initial flat modes propagated into the interaction region must 
be given to a rather high adiabatic order of approximation. 







1 Introduction 



Exact gravitational plane waves are very simple time dependent plane symmetric solutions 
to Einstein's equations Yet, as a consequence of the non linearity of General Relativity 
these solutions show some nontrivial global features, the most conspicuous of which is the 
presence of a non-singular null Cauchy horizon (a Killing- Cauchy horizon, in fact). This 
horizon may be understood as the caustic produced by the focusing of null rays 0. The 
inverse of the focusing time is a measure of the strenght of the wave, thus for an Einstein- 
Maxwell plane wave such inverse time equals the electromagnetic energy per unit surface of 
the wave. This makes exact plane waves very different from their linearized counterparts, 
which have no focusing points and admit a globally hyperbolic spacetime structure. One 
expects that exact plane waves may be relevant for the study of the strong time dependent 
gravitational fields that may be produced in the collision of black holes |3], [| or to represent 
travelling waves on strongly gravitating cosmic strings . In recent years these waves have 
been used in classical general relativity to test some conjectures on the stability of Cauchy 
horizons || [7|, and in string theory to test classical and quantum string behaviour in strong 
gravitational fields || || |10| . Their interest also stems from the fact that plane waves are a 



subclass of exact classical solutions to string theory [ITT], O, ITS 



When plane waves are coupled to quantum fields the effects are rather trivial since they 
produce neither vacuum polarization nor the spontaneous creation of particles, in that sense 
these waves behave very much as electromagnetic or Yang-Mills plane waves in flat spacetime 
fill , H1I- Still the classical focusing of geodesies has a quantum counterpart: when quantum 
particles are present the quantum field stress-energy tensor between scattering states is 
unbounded at the Cauchy horizon, i.e. where classical test particles focus after colliding 
with the plane wave |16). This suggests that the Cauchy horizon of plane waves may be 
unstable under the presence of quantum particles. The classical instability of the null Cauchy 
horizons of plane waves is manifest when non-linear plane symmetric gravitational radiation 
collides with the background wave, i.e. when two plane waves collide. In this case the 
focusing effect of each wave distorts the causal structure of the spacetime near the previous 
null horizons and either a spacelike curvature singularity or a new regular Killing- Cauchy 
horizon is formed. However, it is generally believed that the Killing-Cauchy horizons of the 
colliding plane wave spacetimes are unstable in the sense that "generic" perturbations will 
transform them into spacelike curvature singularities. In fact, this has been proved under 
general plane symmetric perturbations JD]]. Also exact colliding plane wave solutions with 
classical fields are known that have spacelike curvature singularities and which reduce, in 
the vacuum limit, to colliding plane wave solutions with a regular Killing-Cauchy horizon 



Colliding plane wave spacetimes are some of the simplest dynamical spacetimes and, as 
such, they have been used as a test bed for some problems in classical general relativity such 
as the, just mentioned, stability of the Killing-Cauchy horizons, or the cosmic censorship 
hypothesis J7). Note that in a colliding plane wave spacetime with a Killing-Cauchy horizon 
inequivalent extensions can be made through the horizon, and this implies a breakdown 
of predictability since the geometry beyond the horizon is not uniquely determined by the 
initial data posed by the incoming colliding plane waves. The singularities in colliding wave 
spacetimes are also different from the more familiar cosmological and black hole singularities 
which are originated by the collapse of matter since they result from the non linear effects of 
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pure gravity. The type of singularities also differs in the sense that these are all encompassing, 
i.e. all timelike and null geodesies will hit the singularity in the future. 

In a previous paper |T9J which we shall refer to as paper I, we made a first step in the study 
of the interaction of massless scalar quantum fields with a gravitational background which 
represents the head on collision of two linearly polarized shock waves followed by trailling 



gravitational radiation which focus into a Killing- Cauchy horizon ||], 21|. The spacetime 



is divided into four distinct regions: a flat space region (which represents the initial flat 
region before the waves collide), two single plane wave regions (the plane waves before the 
collision) and the interaction region which is bounded by the previous three regions and a 
regular Killing- Cauchy horizon. The interaction region is locally isometric to a region inside 
the event horizon of a Schwarzshild black hole with the Killing- Cauchy horizon corresponding 
to the event horizon. The presence of the Killing horizon made possible the definition of 
a natural preferred "out" vacuum state |22| and it was found that the initial flat vacuum 



state contains a spectrum of "out" particles. In the long-wave lenght limit the spectrum 
is consistent with a thermal spectrum at a temperature which is inversely proportional to 
the focusing time of the plane waves. Of course, the definition of such "out" vacuum is 
not possible when we have a curvature singularity (i.e. in the "generic" case) instead of an 
horizon, whereas a physically meaningful "in" vacuum may be defined in all colliding plane 
wave spacetimes. 

In this paper we compute the expectation value of the stress-energy tensor of the quan- 
tum field near the horizon in the initial flat space vacuum. The expectation value of the 
stress-energy tensor is the relevant observable in the quantization of a field in a gravita- 
tional background since it is the source of gravity (to be written to the right hand side in 
the semiclassical modification of Einstein's equations). These semiclassical equations are in- 
terpreted as dynamical equations for both the quantum field and the gravitational field and 
determine the backreaction of the quantum field on the spacetime geometry. We find, not 
surprisingly, that the stress-energy tensor is unbounded at the horizon. The specific form of 
this divergence suggests that when the backreaction is taken into account the horizon will 
become a spacetime singularity, i.e. the Killing- Cauchy horizon is unstable under vacuum 
polarization. Note that this is a non perturbative effect, it is the result of the nonlinearity 
of gravity, since gravitational waves in the linear approximation do not polarize the vacuum. 
In fact the vacuum stress-energy tensor of a quantum field in a weakly inhomogeneous back- 
ground was computed by Horowitz |23 ] , and it is easy to see that such tensor can be written 
in terms of the linearized Einstein tensor only |24], which vanishes for gravitational waves. 

The non perturbative evaluation of the expectation value of the stress-energy tensor of 
a quantum field in a dynamically evolving spacetime is generally a difficult task. So far 
the most relevant calculation in this respect is, probably, that of the expectation value of 
the stress-energy tensor in the Unruh vacuum state of a blak hole near its horizon [25, |26| . 



Note that although such calculation is done in the extended Schwarzschild spacetime rather 
than in the dynamical spacetime describing the gravitational collapse, the Unruh vacuum 
describes in the extended spacetime the vacuum state with respect to modes which are 
incoming from infinity in the collapsing spacetime. 

Even when the exact modes of the quantum field equation are known it may not be 
possible to perform the mode sums in order to get the quantum field two point function or, 
more precisely, the Hadamard function, which is the key ingredient in the evaluation of the 
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stress-energy tensor. In our colliding wave spacetime we do not even know the exact solution 
of the modes in the interaction region (a similar situation is produced in the Schwarzschild 
case). Fortunatelly the geometry of the colliding spacetime is such that the initial modes 
which come from the flat region are strongly blueshifted in their frequency in the interaction 
region near the horizon. This makes the use of the geometrical optics approximation and, 
more generally, of a systematic WKB expansion of the modes in the interaction region a very 
useful tool. Here we should mention that the blueshift of the initial modes in the interaction 
region is not exclussive of the particular colliding wave spacetime that we take, it is a general 
feature in colliding wave spacetimes and it is due to the focusing produced by the initial 
plane waves. This fact has also been exploited, for instance, to compute the production of 



particles in the so called Bell-Szekeres colliding spacetime |27j . 

The plan of the paper is the following. In section 2 the geometry of the colliding plane 
wave spacetime is briefly reviewed. In section 3 the mode solutions of the scalar field equation 
are given for the four different regions of the spacetime, it is only in the interaction region 
that exact solutions for these modes cannot be found. In section 4 the geometrical optics 
approximation is used to relate the parameters of the modes in the initial flat region with 
the parameters of the modes propagated to the interaction region, in this way a physical 
meaning for these parameters is found. In section 5 the initial flat spacetime modes which 
are propagated in the interaction region are given in terms of a WKB expansion up to 
the adiabatic order four. This is the order which will be required for the regularization 
of the Hadamard function. Then in section 6 the point splitting technique is reviewed 
for the computational purposes of this paper and in the short section 7 the renormalized 
stress-energy tensor is computed in the single plane wave regions, this is not necessary since 
we know from the literature [|15| , |28| that it vanishes, but this calculation is a simple and 
illustrative application of the point splitting regularization technique. In section 8, which 
is the core of the paper, this technique is used to regularize the Hadamard function and 
calculate its value by a mode sum near the Killing- Cauchy horizon. Finally, in section 9 the 
expectation value of the stress-energy tensor near the horizon is calculated. A summary and 
some consequences of our results, such as the backreaction problem, the quantum instability 
of the Killing-Cauchy horizon and the generality of these results are discussed in section 10. 
In order to keep the main body of the paper reasonably clear, many of the technical details 
of the calculations, as well as a short review on the algebra of bitensors, have been left to 
the Appendices. 



2 Geometry of the colliding plane wave spacetime 

Here we recall the main geometric properties of a spacetime describing the head on collision 
of two linearly polarized gravitational plane waves propagating in the ^-direction and forming 
a regular Killing-Cauchy horizon; further details can be found in paper I. This spacetime 
has four regions (see Fig. 1): a flat region (or region IV) at the past, before the arrival of 
the waves, two plane wave regions (regions II and III) and an interaction region (region I) 
where the waves collide and interact nonlinearly. The geometry of these regions is described 
by the following four metrics, which are solutions of the Einstein's field equations in vacuum 
and are written in coordinates adapted to the two commuting Killing vector fields d x and 

dy, 
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,2 a t r r-. • / \l2 , , 1 - Sill (« + f ) , 2 

as T = 4L1L2 1 + sm [u + v)\ dudv ; -ax — 

1 L v n 1 + sin [u + v) 

— [1 + sin (it + u)] 2 cos 2 (it — v)dy 2 , 



a; 



<is 2 n = 4:LiL 2 (1 + sin ufdudv <ix 2 — (1 + sinit) 2 cos 2 it<ii/ 2 , (2) 

1 + sin it 

ds 2 m = ALiL 2 (1 + sini;) 2 dudv ; — dx 2 — (1 + sini>) 2 cos 2 i;<i?/ 2 , (3) 

1 + sm v 

ds 2 w = 4LiL 2 dudv — dx 2 — dy 2 , (4) 

where it and v are two dimensionless null coordinates (v + u is a time coordinate and t> — u a 
space coordinate) and Li, L2 are two arbitrary positive lenght parameters, which represent 
the focusing time (i.e. the inverse of the strength) of the plane waves. The boundaries of 
these four regions are: {u = 0, v < 0} between regions IV and II, {v = 0, u < 0} between 
regions IV and III, {v = 0, < u < tt/2} between regions II and I and {u = 0, < v < 7r/2} 
between regions III and I. At these boundaries the matching of the metrics is such that the 
Ricci tensor vanishes, = (i.e. we have a vacuum solution in the entire spacetime). 

At the boundaries u = tt/2 and v = tt/2 on regions II and III, respectively, the determi- 
nants of the respective metrics vanish, these are the points of focusing of the plane waves 
and are coordinate singularities. The causal structure of these regions is best described with 
the use of appropriate coordinates (harmonic coordinates). In these new coordinates the 
boundaries u = tt/2 and v = tt/2 are seen to be spacetime lines. 

Region I is locally isometric to a region of the Schwarzschild metric bounded by the event 
horizon. This is easily seen with the coordinate transformation, 

t = x, r = M[l + sin(it + v)] , ip=l + y/M, 6 = vr/2 - (it - v ) 
where we have defined M = y/L^L^. Then the metric (|l|) becomes 

ds i = - ' dr 2 -f^-l^j dt 2 - r 2 (d6 2 + sinW) , 

which is the interior of the Schwarzshild black hole. The surface u + v = tt/2 corresponds 
to the black hole event horizon. The boundary v — corresponds to r = M(l + cosd) 
and u = corresponds to r = M(l — cos#), these are the boundaries of the plane waves. 
These boundaries join at r = M (spacetime point of the collision) and also at the surface 
u + v = 7t/2dX8 = and 9 = tt. This region of the black hole interior is outside the 
singularity r = and thus the interaction region has no curvature singularities. The above 
local isometry is not global however, the coordinates 9 and are cyclic in the black hole 
case but in the plane wave case, —00 < y < 00 and —00 < v — u < 00. 

As in the the Schwarzschild case it is convenient to introduce a set of Kruskal-Szekeres 
like coordinates to describe the interaction region, because the (it, v, x, y) coordinates be- 
come singular at the horizon. Since the Klein-Gordon equation can be separated in the 
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new coordinates, these will play an important role in the quantization. First we introduce 
dimensionless time and space coordinates (£, 77) 

£ = u + v , 77 = v — u, (5) 

with the range < £ < 7r/2, — 7r/2 < rj < tt/2. Then we introduce a new time coordinate 
£* related to the dimensionless time coordinate £ by 

£* = 2Mln (^^f] - M (sine - 1) , (6) 
\ 2cos 2 4 J 

and a new set of null coordinates U = £*— x, V = £*+x. Note that the transversal coordinate 
x, which behaves badly at the horizon, appears in this coordinate transformation. Finally, 
we define, 

U' = -2Mexp (-U/AM) < 0, V = -2Mexp (-V/AM) < 0, (7) 
and the metric in the interaction region ([[]) can be written as, 

2p (l-sm£)/2 

dsj = - -dU'dV - M 2 (l + sinO 2 dri 2 - (1 + sinf) WVy 2 , 

(1 + sin 4) 

with 

rr'i/' g »/2 cos 2 e /sin£-l\ ^ / x \ 

U V = 8M -exp , — = exp — — . 

1 + sin^ y \ 2 /' V ^\2MJ 

The curves £ = const, and a; = const, are, respectively, hyperbolae and straight lines through 
the origin of coordinates (U' = V = 0). The Schwarzshild horizon (which is a Killing- Cauchy 
horizon for our spacetime) corresponds to the limit of the hyperbolae when £ — ► 7r/2 i.e. 
V = or U' = 0. Notice that the problem with the transversal coordinate x at the horizon 
is that all the lines x = const, go through the origin of the (U 1 , V) coordinates, so that all 
the range of x collapses into the point V — U' — 0, whereas the lines U' — and V — 
represent x = —00 and x = 00 respectively. One should recall that we have not represented 
the coordinate x in Fig. 1 where only the (u, v) coordinates are shown, a; is a transversal 
coordinate perpendicular to the propagation of the waves. 

To understand the global geometry of the spacetime a tridimensional picture like Fig. 
3 of paper I, is helpful. It represents the boundary surfaces between the different regions 
in terms of the appropriate nonsingular coordinates adapted to each region. The spacetime 
lines £ (u = tt/2) and CJ (v = tt/2) are identified with points V and V respectively, these 
points are known as folding singularities, and are avoided by all null geodesies except for a 
set of null measure. There is no spacetime beyond the lines £ and £' (see Fig. 1). 

3 Mode solutions of the field equation 

In order to quantize a field in our background spacetime and to compute the field two-point 
function which will be needed later, we need the mode solutions of the field equations. Here 
we consider a massless scalar field 0, the Klein-Gordon equation is 
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□0 = 0, (8) 

where □</> = (— g)~ 1 ^ 2 [(— 9) X ^^A,^ an< ^ 9 * s t ne determinant of the metric g^ u . Note that 
since the curvature scalar R is zero in our case, the above equation is valid whatever coupling 
to the curvature we may consider. The mode solutions of these equations are different in the 
four spacetime regions we have. In paper I details of the equations in the different regions 
and their solution by separation of variables were given. Here we summarize the main results 
and introduce some new ones which will be of use later. 

Our vacuum state is the "in" vacuum, this is the physically unambiguous vacuum defined 
in the flat region IV before the plane waves arrive. In this region a complete set of positive 
frequency modes with respect to the timelike Killing vector d u / +v i is given by: 

uf(u, V, X, y) = 1 e -i2L 2 k-V-i2L lk+ u+ik x x+ik y y ^ (g) 

^2k_(2nf 

where u and v are the two dimensionless null coordinates related to a physical null coordi- 
nates u' and v' by u' = 2L\u, v' = 2L2V. The labels k x , k y and k- are independent separation 
constants for equation (|8]) and the label k + is determined by the relation, 

Ak + k_ = k 2 x + k 2 y . (10) 

It was shown in paper I that these modes are well normalized on the hypersurface {u = 
0, v < 0} U {t> = 0, u < 0}. The labels k x and fc_ are continuous but k y is discrete if we take 
a cyclic spacetime in the y-direction (this is convenient but not necessary), then we identify 
k y with mM~ x where m is an integer. 

These "in" modes, u™, which define our vacuum state, have to be propagated through 
the spacetime. Once we know the "in" modes in a certain region we use their expression on 
the boundaries with the next region as boundary conditions for the "in" modes in this next 
region. These will be given in terms of a complete set of solutions of the field equation in 
such region. In the plane wave regions II and III the "in" modes are easy to find: 



U%(x) = 1 e ik x x+ikyy 

/2kj2n) 3 



L . e - l 2L 2 / £ _,- J A(«)/(2L 2fc _) ; RegionII 



cosu 

1 



111 



e -i2L 1 k + u-iA(v)/(2L 1 k + )^ Region m 



COSf 

where A(x) = kff(x) + k\ g{x) and where the two fuctions f(x) and g(x) are 

„, . (1 + sinx) 2 .15 15 . . . . 

j(x) = (9 — smxH cosx x — 12, g(x)=tanx. (12) 

2 cos a; 2 2 

We use for simplicity the dimensionless labels k\ = Mk x and k 2 = Mk y where k x , k y and fc_ 

have the same meaning as in the flat region IV because these expressions for the solutions 

of the Klein-Gordon equation in region II and III match smoothly (i.e in a continuous and 

differentiable way) with the respective solutions @ on the boundary between regions II and 

IV, i.e {u = 0, v < 0} and the boundary between regions III and IV, i.e {v = 0, u < 0}. 

In the interaction region (region I), however, this kind of mode propagation is not that 

simple. First, we note that equation (|8D in this region can be separated as 
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<KZ,V,x,y) = e ik * x+ik » y i> ak Mv*k y (v), (13) 
where coordinates £, r\ are related to the usual null coordinates u, v, by (|^) and the differential 
equations for i[>ak x (0 an d fak y (v) are (dropping the indices) 

V>« - (tanO^ + fa + ^ (l + smO \ 
\ cos^4 ) 

tp m - (tan 77)^ + ( a — h\ — —J <p = 0, (15) 

V COS z 1] J 

where a is a dimensionless separation constant and k x , k y (or k\ = Mk x , k 2 = Mk y ) are the 
same labels as in regions IV, II or III. 

The solutions of fll5[) can be expressed in terms of associated Legendre functions and in 
the cyclic case they can be written in terms of spherical harmonics, Y l m (y/M, 77/2 — 77), with 
discrete labels related to a and k 2 simply by a — 1(1 + 1) and k 2 = m. The solutions of fll4j) , 
however, cannot be written in terms of known functions. If we use the variable instead of 
£ defined in (^), i.e. d£ = cos£ M _1 (l + sin£)~ 2 and introduce a new function 7 defined 
by 27 = (1 + sin£) ip, equation ([14]) becomes 

7,^*+"i7 = °> (16) 

where 



ul(0 = M~ 



cos 2 £ 2 
kl+ M : -A-u [1(1 + 1) + 



1 + sin0 4 



1 + sin £ 



M~ 



k\ + V 1 (0 



(17) 



Note that a new function Vi(£) has been defined. Equation flilf ) can be solved by a WKB 
method and the solutions 7(£*) can be expanded up to any adiabatic order as 



7(D 



w x (0 



where Wi(£) is given, up to adiabatic order 4 (this is the order we need) 



^1(0 = ^ + — + — + — + — + -9 + — > 



Wi CUV CUT CUi 



where, 



by 



(19) 



8 M 2 ' 



Bo = — 



b_vi_ 

35 M 4 



(20) 



J_ Vj_ 1 281^ + 19 If _ 221 Vg_t^ 1105 ^ 4 

4 ~32M 2 ' 4 ~ 128 M 4 ' ~~ 258 M 6 ' 4 ~ 2048 M 8 ' 

an overdot means derivative with respect to £*, and A n , S n , ... denote the n adiabatic terms 
in Wi(£). From the two exponential factor signs in the WKB solution (|i~8|) , we choose the 
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negative sign for consistency with the boundary conditions imposed by ([□]), as we will see 
later. This gives positive frequency solutions with respect to the timelike vector d/d^* on the 
horizon since these solutions are asymptotically proportional to exp(— i\k x \^*). A complete 
set of solutions of the Klein-Gordon equation in the interaction region is thus given, in the 
cyclic case, by 

Mx) = ]° jL -e^Jo^'^Yriy/M^^-r)), (21) 



where C is a normalization constant which can be easily calculated imposing that these 
solutions are well normalized on the horizon (£ = 7r/2). We will find its value in a next 
section. 

Finally, the "in" modes in the interaction region, with the boundary conditions imposed 
by (PD) can be written as, 

7j, x,y) = Y, Ci <M£, V, x, y) , (22) 
i 

where the coefficients Ci depend on I and on the separation constants used to label the modes 
in region II, i.e k x , k y , k_. We devote the next two sections to find an explicit expression 
for 



4 Geometrical optics 

The solutions of the Klein Gordon equation in the flat and single plane wave regions, (|9|) and 
(|Tl|) respectively, can be understood in terms of geometrical optics because the geometrical 
optics approximation is exact in these regions, i.e. the orthogonal curves to the surfaces 
of constant phase (the rays) follow exactly null geodesies. This will allow to "localize" the 
Klein-Gordon solutions in the spacetime and to determine which solutions will be relevant 
in later calculations. The geometrical optics also gives information on the meaning of the 



parameters which label the "in" modes (JTTJ) , i.e. k±, k x , k y , and the parameters k x , m, I, 
which label the modes (^Tj) in the interaction region, and on their connections. 

In the flat region the normalized "in" modes ([J) are purely flat modes which have positive 
frequency with respect to the timelike Killing vector d t = d u r +v r. The energy of these modes, 
Eiy, may be defined as the eigenvalue of d t , d t u^ = —iEjyu^ . In terms of the null 
momenta k± it is, 



E w = k + + k„. (23) 

We can now proceed in the same way in the single plane wave regions. In region II, the 
normalized "in" modes take the form flHD , but now, instead of a timelike vector, dt, we 
have a null Killing vector d v >, and ( |TT]) are its eigenfunctions with positive eigenvalue k_, 
i.e. d v i = —ik-Ufr . Similarly, in region III, k + are eigenvalues of d u >. In this 

case, however, these eingenvalues cannot be directly interpreted as energies. The mode 
solutions in the plane wave regions II and III have the form C e lS , where S = constant 
define the surfaces of constant phase of the modes and it satisfies exactly the null condition 
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g^S^S^ = 0. Thus, we might define the energy of these modes as the variation of S in the 
direction of the vector —d t , i.e. £?n/m = — d t Sii/m. From (|TTD , we identify 



S n = -2L 2 k_v — — h k x x + k y y, 

2L 2 k_ 

and Sjjj, by a similar expression. The tangent vector fields to the congruencies of null 
geodesies, expresed in (u, v, x, y) components are, in region II, 

and a similar expression for Vj H in region III. 

By integration we obtain a congruence of null geodesies parametrized by their momenta 
k±, k x , k y and their initial positions v , x , y at u = in region II and Uq, x , y at v = 
in region III. Each congruence at fixed values of the momenta represents the set of rays 
orthogonal to the surfaces Sn/m = constant, so that these rays "localize" in the spacetime 
the mode labeled by k±, k x , k y . A simple inspection of the expression for the tangent vector 
fields above shows that in the plane wave region II the rays for each mode are peaked near 
the line C(u = vr/2) and reach the interaction region at points close to the folding singularity 
V (u — 7r/2, v — 0) if the momenta satify the conditions, 

2L 2 k_ > 1, 16A^ + Af < 2L 2 k_, (24) 
where the second condition comes from the property, 

f(u)kl + g[u)k\ < (16Jfe? + k\) cos~ 2 m, 



of the functions / and g, defined in ([12]). Similar conditions hold for modes in the plane 
wave region III, by the simple substitution of 2L 2 k_ by 2L\k + , i.e. the rays for each mode 
are peaked near the coordinate singularity line £' (v = tt/2) and reach the interaction region 
at points close to the folding singularity V' (u = 0, v = 7r/2) if 

2Lik + > 1, 16kl + kj < 2L x k + . (25) 
Due to the relation (10), the two conditions (^3j) and (25) are mutually exclusive. These 



geometrical properties of the "in" modes indicate that any calculation involving "in" modes 
in regions close to the folding singularity V, requires only modes in region II with momenta 
satisfying (El) and no modes in region III. Equivalently any calculations in regions close to 
the folding singularity V, requires only modes in region III with momenta satisfying (|25|) 
and no modes in region II. Note that the larger is the energy Ejy of the "in" modes in the 
flat region, the closer their rays get to the lines £ or £', and this means that the "in" modes 
are blueshifted towards these lines and consequently towards the horizon. In fact, this is a 



general property of plane waves |T6| and it will become very important in the evaluation of 
the vacuum expectation value of the stress-energy tensor when summations over "in" modes 
are performed. 

Going back to the energy of the "in" modes in the plane wave region II we have, according 
to the previous definition, En = —d t Su = + (f(u)kl + g{u)ky) (4A;_) _1 , which is not 
constant, because dt is not a Killing vector, and cannot be interpreted properly as the energy 
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of the mode. However, close to the origin u = v = 0, i.e. just where the plain waves collide, 
En coincides with Ejy, since 



En\u,v~o = k- + + °(« 2 ) = k - + K + 0(u 2 ) ~ E IV , 

and close to the line C we obtain, in powers of cosw, 

E u = k_+ 16fc f, + ky cos" 2 it + 0(cos° u). 

Ak_ y ' 

However, if we recall that for the "in" modes close to the line L in region II the inequalities 
([24]) are satisfied, then there is a large interval of the variable u for which, 

* , 2 cos~ 2 m< 1. 26 

In this interval Eu — k- — Eiy, which is the eigenvalue of the null Killing vector d v r, thus 
it can be interpreted unambiguously as the energy of the "in" modes in region II. Similarly 
for "in" modes close to de line CJ in region III there is a large interval of the variable v for 
which, 

Meos-. t ,«l. (27) 

and their energy is unambiguously defined by Em ~ k + ~ Eiy. We will see in the next 
section that these intervals, fl2"6| ) and (E7|), play an important role. 

In region I, the mode solutions of the Klein-Gordon equation are given by fl2ip . Although 
the geometrical optics approximation is not exact in this region we can still make good 
use of it, since this approximation is compatible with the WKB expansion of equation 
(|i6|), which depends on the dimensionless timelike coordinate £ = u + v (or £*). The 
energy of these modes can be defined by the variation of the surfaces of constant phase, 
S = -J$ Wi (£')*;'* + ik x x + ik y y, along the timelike vector M~ l i.e, E T = —M~ l d i S ± . 
If we consider the range of values p = |fci|L(cos£) _1 <C 1 in analogy with ( p6|) and ( p7[) , 
Ei ~ M _1 [l + 0{p 2 )} L~ l . In other words, the parameter M^L" 1 = M~ l [l(l + 1)] 1/2 , that 
labels the mode solutions of the Klein-Gordon equation in the interaction region, takes the 
role of an energy. Recall that this energy is positive just because we have choosen the minus 
sign in the WKB solution (fT8p. This suggests that we can take the value for the "in" modes 
in the interaction region close to the horizon, w™, as a single function 0^ ( pi]) instead of the 
linear superposition (p2|) , provided that the energy of fa, i.e Ei equals to En and to Em, at 
the boundaries of regions II and III with region I, respectively. In other words we can take 
it™ ~ 0^ where 4>k is labeled by a parameter L satisfying, 

Ei ~ M^L' 1 oc + k+. (28) 

In the next section we will see that this identification is really possible and we will find 
the proportionality constant. Note that the "in" modes which are relevant at the boundary 
between region II and region I satisfy 2L 2 k_ ^> 1 and 2Lik + ~ 0, so that Ei ~ Eu and 
L oc L 2 k-, and the modes which are relevant at the boundary between region III and region 
I satisfy 2Lik + ^> 1 and 2L 2 k^. ~ 0, then Ei ~ Em and L oc Lik + . 
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5 Adiabatic expansion of the "in" modes 

We seek a solution of the Klein-Gordon equation near the horizon that satisfies the boundary 
conditions imposed by the "in" modes (|TI|), which have been propagated from the flat region 
through regions II and III. The problem is formally solved using a Bogoliubov transformation 
provided {u k n }, i.e. the "in" modes, and the set of modes fl2lD , {4>k'(%)}, of the field equation 
in the interaction region are complete and orthonormal, so that we can write, 

u k ( x ) = J2 a kk> [4>k'(x) + (3 kk r(f)* kl (x)] . 
k' 

All the information on the boundary conditions is contained in the Bogoliubov coefficients, 



akk' = (u k n (x), 0fc/(a;))s, (3 kk > = -(u k n (x), (29) 

where the inner product is evaluated on the boundaries of regions II and III with region I 
and it is given, see paper I, by, 



'1, V2 



) = — i / dxdy 



tt/2 



COS 2 M (0i<9 u 02 



/"7r/2 «-> 

du+ cos 2 v (4>id v <f)2) 

v=0 JO 



dv 



u=0 



(30) 



The coefficients a^} and /3j£) can be calculated using (p9|) , (|30| ) and the expressions 
(PI) for the functions <fi k >(x) and u k n , respectively. The superscript (A) indicates that the 
values of <p k >{x) are known up to a certain adiabatic order A. We have, 



41' = 1 



Mv^F (-l) m 5(k x - k' x ) 5 mm , A$, = -iM^5(k x + k' x ) 5 m( _ m0 fiW 



(A) 



where ^4^) and B^l are given by. 



(A) 



•^kk' 



2C r/2 



\fkl Jo 
2C r/2 



JO 



G?W COS K 



<iw cos f 



4j 4( U )/(2L 2 fc_) „ / j jj* Wi(€')d£" 

l + sin£ ^(f) 



COSM 



<9 U 



+ 



l + sin£ 



COS 1> 



d v 



v=0 



u=0 



(31) 



»(A) 
°kk' 



2C r/ 2 



\fkl Jo 

2C r/2 



du COS K 



dv cos u 



iA(u)/(2Lafc_) „ ^ e ~</ o e Wi(Od£' 



COSM 



<9 U 



1 + sin £ 



+ 



<A(«)/(2L 1 fe f ) „ ^) e -i/ o e Wi(€'K'* 



where <p(rj) is a shortened notation for 
<p( V ) = YT(P,*/2-Ti) 



cosv 



l + sin£ Jwjg) 



v=0 



u=0 



'21 + 1 (i - m y: 

47T (l + m)\ 



1/2 



^(-sinry), 



(32) 



(33) 
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and where Wi(£) is known up to a certain adiabatic order A. 

The problem now is to evaluate these integrals. Since all of them are of the type, 
/ dx /(x)exp we search for the intervals of the integration variable where the modulus 

of the integration function, f(x), is much larger than the variation of its phase 8(x), i.e 
6(x) -C f(x). In particular, we search for the phase stationary points, i.e. the values of the 
integration variable for which 6(x) = 0. 

Using that £ = u + v, rj = v — u and defining star variables x* by dx* = dxM(l + 
sin a;) 2 / cosx, the phases of the integration functions in ( pT|) and (|32| ) can be written, re- 
spectively, as 



A(u) 



± 



W 1 (u')du'*, 



A(v) 



± / Wi(v')dv' 
J o 



(34) 



2L 2 k_ Jo v ' 1 2L x k 
We can study the problem qualitatively using the value of W\{£) up to adiabatic order zero 



(o) 



1 

M 



k\ + 



cos 2 u 



1 



+ 



l + sinu) 4 Vl/ 2 1 + sinu 



1/2 



(35) 



where we have used the notation 1(1+1) = L~ 2 . We have already seen in the previous section, 
using geometrical optics, that the relevant modes involved in any calculation near the horizon 
(£ = 7r/2) are those which satisfy (|26"D in region II (or (|27 ) in region III) and this implies (in 
region II) that the labels L^ 1 and L2&-, which parametrize the Klein-Gordon modes (plf) 
in the interaction region, and the "in" modes ( ITT] ) are related. Let us assume that they are 
proportional, then it is easy to see that expressions fl34|) have stationary points only when 
the minus sign is involved and these stationary points are such that L~ 2 cos 2 u is of the order 
of 16/c 2 + k\. This means, in particular, that L _2 cos 2 u is of the order of k\ = m 2 ; but for 
these values, the associated Legendre function -P/™(— sinw) has an oscillating behaviour that 
we must take into account in order to evaluate the phase stationary points. Note that for the 
range of values 16k 2 + k 2 ^> L~ 2 cos 2 u, the associated Legendre polinomial -P; m (— sin-u), with 
k 2 = m, behaves as cos m w and in that case we have integrals / dx f(x)exp [i6(x)] without 

~ cos~ 2 u, which lead to negligible 



,m-l 



u -C 



[u 



stationary points and with f(u) ~ cos" 
contributions. 

In the following we will find the expressions of the Klein-Gordon solutions <pk in the 
interaction region for the range of values such that |/ci|L(cos^) _1 1 and \k 2 \L(cos r])^ 1 <^ 
1. We have seen that solutions <pk can be separated in two functions ip(£) and (p(r]), i.e. 
( fik(£,,v) = ^(Ov^l 7 ?)) which satisfy the differential equations (0) and (|T5[). The first of these 
equations can be adiabatically solved up to any order, see (|I8|), and the second one, although 
exactly solvable in terms of associated Legendre functions, can also be WKB solved up to 
any adiabatic order in the range of values such that | A:2 |i^(cos ^y) - 1 1. For this we define 
a new coordinate 77* by drj* = drj M(cos^) _1 in (|15|) and obtain, 



where 



1 

M 2 



ip tV * v * + uj 2 if 
-kl + 



2 , COS 2 Tj 



L 2 



= 0, 
1 

M 2 



-k 2 2 + V 2 ( V ) 



(36) 
(37) 



which is WKB solved for the values |A; 2 |-^(cos 77) 1 <C 1 as 
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Co 



W 2 (v) 



j / " www* 



where W 2 (rj) can be written up to adiabatic order 4 in exactly the same way as (19|), 



(38) 



W 2 (r]) 



A 2 

^ + -3 + 

U> 2 



5 4 



a, 



-Dzi 



/ >9 
W 2 



11 ' 



(39) 



where the expression of A 2 , B 2 , ■ ■ ■ are formaly the same expressions ( |20| ) for A 2 , B 2 , • • -, but 
with Vi(£) substituted by V^??) (defined in (|37|)), and with the overdot meaning derivative 
with respect to rj*. The two arbitrary constants C\ and C2 will be fixed imposing that the 
solution (|38|) is asymptotically equivalent to the associated Legendre function P/™(— sinr/) 
with L 2 = [1(1 + 1)] _1 and m = k 2 , for the values Ik^L^cosrj)^ 1 <C 1 (see Appendix A). 

Before going through, we need to introduce some useful definitions to relate the adiabatic 
expansions of Wi(£) and W 2 (rj). Let us define, 



q 2 = (16kf + k\ 



V 1 



\§k\ q 2 , 



U 2 



,2 2 
" 2 H 1 



note that in practice, and without loss of generality, we can take q 2 < 1, because the case 
q 2 ^> 1, which corresponds to k\ ~ k 2 ~ 0, is equivalent to consider L <C 1 in the equations 
(H) and (H). Note also that < U, V < 1, and that by constuction, V 2 + U 2 = 1. With 
these definitions we can write, 



M 2 u;j 



V 2 



cos 2 £ 




M 2 u; 2 



U 2 

— + 



cos 2 r/ 



16 g 2 (l + sin£) 2 V L2 ' 1 + sin ^ / 

For parameters that satisfy IfcilL(cos^)^ 1 <C 1 and Ik^L^cosr])^ 1 <C 1, we introduce new 
functions p(£) = Lg _1 (cos^) _1 and ^1(77) = Lq~ l (cos 7y) _1 , so that p <C 1, pi <C 1. 

We can now expand (^1|) and ([□]) in powers of p and pi and substitute into ( |3T] ) and 
(0). Then the dominant integrals connecting the "in" modes between regions II and I can 
be written as (see Appendix A for details). 



2M / duexp 



du 



' (V 2 U r 

p2 U (1 + sinM)4 + T 



L 



Uk_ 



0( P 4 



L)+0(p 4 



5 (L 2 k- 



2M 



L~ 



(40) 



The dominant integral connecting the "in" modes between regions III and I is obtained from 
this by changing L 2 k_ by L\k + . This means that only the a.f.y connection coefficients (|29| ) 
are relevant in the propagation of the "in" modes through the boundaries between regions 
II and III with region I. In other words, the set of modes of the Klein-Gordon equation 
for region I correspond to the same initial vacuum state for calculations involving points 
near the horizon £ = tt/2. This important result, as we have pointed out in the geometrical 
optics analysis, can be seen as a consequence of the blueshift that the "in" modes suffer by 
a single plane wave. Although /c_ (or k + ) is a continuous label and L is a discrete label, this 
should not be a problem since we can use the wave packet formalism to discretize k_ []3~0[]. 
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6 "Point splitting" regularization technique 



In this section we briefly review, for the computational purposes of the following sections, the 
"point-splitting" regularization technique to calculate the expectation value of the stress- 
energy tensor of a scalar quantum field in some physical state. The stress-energy tensor 
of the field may be derived by functional derivation of the action for the scalar field with 



respect to the metric. When the field is massless and the Ricci tensor is zero it is [51 1 



T, 



[IV 



;i-2£)> 



■9tu>9 



';a(P;P 



2£. 



+ 2^0V a V a 0, (41) 



where £ is the coupling parameter of the field to the curvature. 

To quantize, the field is promoted into a field operator acting over a given Hilbert 
space |3l], |3~3l , <p(x) = J2k a kUk(x) + a\u* k (x), where the standard creation and 



annihilation operators and {uk(x)} is a complete and orthornormal set of solutions of the 
Klein-Gordon equation (|8|). Mathematically the field operator <j>(x) is a point distribution, 
therefore, the quantum version of the stress-energy tensor ( ^Tj ) is mathematically patholog- 
ical because it is quadratic in the field and its derivatives. One possible way to give sense 
to that expression is to note that the formula ( f41~|) can be formally recovered as, 



(T fa/ )(x) = lim V^G^{x,x'), 



(42) 



where G^(x, x') is a Green function of the field equation defined as the vacuum expectation 
value of the anticommutator of the field, and called the Hadamard function, 



G^(x,x') = ({0(x),<Kx')}> =Y,{ u k{xX{x')+u k {x')ui{x)}. 



(43) 



As a product of distributions at different points this is mathematically well defined. The 
differential operator is given in our case by, 



V 



ij,i> 



(1-20 I (V^ V„ + V^) + (2 £ - g, u 1 (V a , V a + V Q V a ') - 

£ \ (V^V, + V^V,0 +tg» v \ (V Q V Q + V Q ,V Q ') . (44) 



However, the above differential operation and its limit has no immediate covariant meaning 
because G^'(x, x') is not an ordinary function but a biscalar and the differential operator 
T)^ is nonlocal, thus we need to deal with the nonlocal formalism of bitensors. In Appendix 
B we give a summary of the main properties of bitensors, how to parallel transport from 
points x to x' and how to average after the coincidence limit, x —>■ x', is taken. 

The above procedure still leads to a divergent quantity since we know that even in flat 
spacetime G^(x, x') has a short- distance singularity and that a "vacuum" substraction has 
to be performed to G^'(x, x') in order to obtain a regularized value. To regularize we assume 
that G {l) (x,x' ), has a short-distance singular structure given by 



S(x, x') 



(4tt) 



;A 1/2 (X,X') 



a(x, x') 



— + v(x, x') lncr(x, x') + w(x, x') 



(45) 
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where a(x, x') is the geodetic biscalar (see Appendix B), A is the Van Vleck-Morette determi- 
nant ||35|| , which is singularity free in the coincidence limit, and where v(x, x') and w(x, x') are 
biscalars with a well-defined coincidence limit for which we assume the following covariant 
expansions, 

oo oo 

v(x, x') — y^ y Vi(x, x')a (x, x'), w(x, x') = ^2wi(x, x')a l (x, x'). (46) 

1=0 1=0 

A Green function expressed in this form is usually called an elementary Hadamard solution, 
the name of which comes from the work of Hadamard on the singular structure for elliptic and 
hyperbolic second order differential equations. Note, however, that this Hadamard singular 
structure is not a general feature of any Green function of the Klein-Gordon equation. In 
other words, although for an extensive range of spacetime and vacuum states, the vacuum 
expectation value of the anticommutator of the field, G^(x, x'), has this singular form, 
this is not a general property. However a theorem states that if G^'(x, x') has the singular 
structure of an elementary Hadamard solution in a neighbourhood of a Cauchy surface of 
an arbitrary hyperbolic spacetime, then it has this structure everywhere |22|, j32fl . As a 
corollary of this theorem, G^(x, x') has this singular structure if the spacetime is flat to 
the past of a spacetime Cauchy surface, as is the case of our colliding plane wave spacetime. 



This and other considerations lead to a proposal by Wald [S3] that any physically reasonable 
quantum state must be a Hadamard state, that is to say, a state for which (x, x') takes 
the short- distance singular structure of an elementary Hadamard solution. 

The coefficients Vi(x,x') and Wi(x,x') can be directly obtained by substitution in the 
differential equation, O x S(x, x') = 0. Recursion relations for Vi(x,x') and Wi(x,x') are 



then obtained [j35] . These relations uniquely determine all the Vi(x, x') coefficients but the 



coefficients Wi(x,x') can be written in terms of an arbitrary term w (x 1 x'). Up to order a, 
v (x, x 1 ) is given by 

v (x, x') = —ai(x, x') — - a 2 (x, x') a H , (47) 



where ai(x, x') and 02 (x, x') are the Schwinger-DeWitt coefficients [jHJ, which in our case 
{R^u = 0) reduce to 



ai (x, x') = --i- R al3 \ R aMv ^u v + 0(a 3 ), a 2 (x, x') = R a ^ s R aPjS + 0(a). (48) 

loll loU 

Only the coefficients Vi(x, x') are related to the singular structure of G^'(x, x') in the coin- 
cidence limit, and they are uniquely determined by the spacetime geometry. This means that 
given any two Hadamard elementary solutions in a certain spacetime geometry, both have the 
same singularity structure in the coincidence limit; therefore given two vacuum Hadamard 
states, |0) and |0), G^(x,x') = (O|{0(x), 0(x')}|O) and G (1 \x, x') = (O|{0(x), 0(x')}|O), 
they have the same singular structure. Their finite parts, however, may differ because the 
two vacuum states are related to different boundary conditions, which are global spacetime 
features. Mathematically this comes from the fact that the term wo(x,x') in the elemen- 
tary Hadamard solution is totally arbitrary, fixing Wq(x, x') we fix a particular boundary 
condition. This suggests a possible renormalization procedure |4^, |3JJ: we can eliminate 
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the non-physical divergences of any :, x') without alterations in the particular physical 

boundary conditions by subtracting an elementary Hadamard solution with the particular 
value w (x,x') = 0, which is the value that corresponds to the flat space case. In other 
words, we define the following regularized biscalar, 

G%\x,x') = G^\x,x') - S(x,x')\ m=0 . (49) 
Then by means of Gg> (x,x') we can construct a (T^J) by differentiation with the nonlocal 



operator (g4J). 

This regularization procedure, however, fails to give a covariantly conserved stress-energy 
tensor since it can be seen that for a massless conformal scalar field (i.e. £ = 1/6), 



V(T*> = - lim V„ (n x , + -R(x')) G B (x,x') = — V M a 2 (x), (50) 



1 

-a«" - a* V" x ' 6" V ~V " y 64^ 2 

where 02(0:) is the coincidence limit of the Schwinger-DeWitt coefficient a 2 (x, x') in (38). 



In particular for a spacetime with null curvature R = 0, such as our colliding spacetime, 
this is true for any coupling. Thus to ensure covariant conservation, we must introduce an 
additional prescription: 



a 2 (x) 



(T^W) = (TZW) ~ 9r- ( 51 ) 

Note that this last term is responsible for the trace anomaly in the conformal coupling, 
because even though (T^ u (x)) has null trace when £ = 1/6, the trace of {T^ u (x)) is given by 
(Tp = -a 2 (x)(16n 

The regularization prescription just given satisfies the well known four Wald's axioms 



2\-l 



37] , |3q , [391 , p3| , a set of properties that any physically reasonable expectation value of the 
stress-energy tensor of a quantum field should satisfy. There is still an ambiguity in this 
prescription since two independent conserved local curvature terms, which are quadratic in 
the curvature, can be added to this stress-energy tensor. This two-parameter ambiguity, 
however, cannot be resolved within the limits of the semiclassical theory, it may be resolved 



in a complete quantum theory of gravity [33] . Note, however, that in some sense this 



ambiguity does not affect the knowledge of the matter distribution because a tensor of this 
kind belongs properly to the left hand side of Einstein equations, i.e. to the geometry rather 
than to the matter distribution. 



7 Vacuum expectation value of the stress-energy ten- 
sor for a single plane wave 

A simple nontrivial example that nicely ilustrates the point splitting technique for the regu- 
larization of the vacuum expectation value of the stress-energy tensor is the particular case 
of a single gravitational plane wave, i.e. the case of regions II and III, (T^u and {T^ u )jjj. 
The result is not new, it is known that the expectation value of the stress-energy tensor in 
the "in" state, i.e. the vacuum state at the flat region before the passing of the wave, is ex- 
actly zero 28]. Let us start with the construction of the Hadamard function G^'(x,x'), 



([43]), if the vacuum state is the "in" vacuum state, we have 
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dk x e ikx ( x - x ')- l2L i [/(«)-/(«')] Z* 00 ^ e **« (y-2/')-i2ii [9 («)-<?(«')] + c . c . ; 



where it/- are the "in" modes (|IT1) of region II. It is easy to compute analytically the above 
"mode sum" since only Gaussian integrals are involved, the result is 



sect* sectt' 



4M 2 (2tt) 2 y/\f(u)-f(u')]\g(u)-g(u')] 

{^M^ + -Mn)T- (52) 

We can now proceed with the point splitting formalism by setting the points x and x' at 
the endpoints of a non null geodesic parametrized by its proper distance. Let us denote by 
f the midpoint on the geodesic, which lies at equal proper distance from x and x', i.e, 

xf = af t {f + e), x? = xP{t), x'» = x^if - e), 

where one should be aware that the notation x, x' may be a little ambiguous because it 
makes reference to the points (u,v,x,y), (u', v', x', y') and also to their third components. 

Now we can easily solve the geodesic equations, since the momenta p x , p y and p v are 
three constants of the motion. The geodesies are a family of curves parametrized by p x , p y 
and p v , 

x = x(u) = x(0) -2M 2 — /(«), y = y(u) = y(0)-2M 2 ^g(u), (53) 

Pv Pv 

,2 M 2 

Pv ' Pv ' ' Pi 

where f(u) and g(u), are given in ([T2p. If we substitute these curves into fl52p using that 
x' = x(u'), y' = y(u'), v' = v(u') and take u, u, u' as u = u(f + e), u — u(f), v! = u(f — e), 
then we only need to expand fl52j) in powers of e up to order e 2 . The calculation is simple 
but tedious and finally one obtains, 

cs i )( rr >\ = LI En z (^) 

K ' >~ An 2 a 160M 8 tt 2 (l + sinw) 10 ' { ' 

where a = £s 2 /2 = 2 e 2 S and s is the proper distance along the geodesic between x and x', 
see Appendix B. 

The regularization follows from the subtraction of the elementary Hadamard solution 

(H) 

1 1 

S(x, x)\ 



r? v 2 M 2 r u 

V = V ( U ) = V (Q) + M 2 ^f{u) + M 2 ^g{u) + -^- / ^'(l + sinw) 2 , (54) 

pi p:, p:, Jo 



'«0 4?r 2 a 
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since the biscalars a 2 (x, x') and ai(x, x') are zero up to order a. Note that, in particular, 
there is no trace anomaly because a,2(x) = 0, thus after regularization there remains a term 
proportional to 

P^ = 4^ = 4^ + 0(e), 
a a 

where we have used that cr M = 2ep^ being p^ a tangent vector to the geodesic connecting x 
with x' such that p^p^ = S. Under application of the differential operator ( [44]) (details are 
given in section 9) one easily gets an expression which involves factor terms of the kind p v p v , 
p v p v PvPfi and PvPvPvPvPftPu- These terms are path dependent in the limit of coincidence 
x — > x' therefore we must introduce the elementary averaging procedure of Appendix B, 
which gives, 

(p v p v ) oc g vv = 0, (pvPvPvPn) oc g vv = 0, (PvPvPvPvP^Pu) oc g vv = 0. 

Thus we recover from P^ ) the known result that the vacuum expectation value of the 
stress-energy tensor of a quantum scalar field, in a single plane wave region is exactly zero, 
(T/j,u)ii = 0, i.e. exact gravitational plane waves do not polarize the vacuum. 



8 Hadamard function in the interaction region 

Here we calculate the Hadamard function G^'(x,x') in the "in" vacuum state in the inter- 
action region, near the horizon. As we have seen in section 5, due to the blueshift of the 
"in" modes towards the horizon, the u™ with large energy values (which are the relevant 
modes in calculations near the horizon) determine the same vacuum state as the modes <fik, 
. This means that near the horizon the Hadamard function in this vaccum state can be 
written as, 

GW(^')^{<(^rM} +cc = E #00} + c - c - (56) 

k k 

Since we do not have an exact expression for the solutions <pk we need to work up to a certain 
adiabatic approximation. This means that we have the inherent ambiguity of where to cut 
the adiabatic series. Fortunately the boundary conditions of the problem impose that the 
asymptotic value for the "in" modes on the horizon is given by, exp (—i\k x \£*), see paper I 
for details. This is just the form for the "in" modes that we used in paper I to calculate the 
production of particles, and is equivalent to cut the adiabatic series at order zero. Although 
for the particle production problem it was enough to cut at order zero, it is not sufficient for 
the calculation of the vacuum expectation value of the stress-energy tensor. This is because 
G^ at order zero does not reproduce the short-distance singular structure of a Hadamard 
elementary solution (fj5 ) in the coincidence limit x — *> x'. The smallest adiabatic order which 



we need, to recover the singular structure of is order four. 

In the mode sum of (KB) we use the shortened notation J2k = I^L dk x Z)£o E^=-d an d 



we note from ([H]) that the sum over m is straightforward using the property of the spherical 
harmonics, 
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m=l 21 + I 

£ YT(7r/2 - 77, y/M)Yr (ir/2 - rf, y'/M) = fl(oose), 

m=— i 

where P/(x) is a Legendre polynomial of order / and G is the angle between the points 
(9, if) and (#', if') on the sphere, i.e, cos0 = cos 9 cos 9' + sin# sin#' cos {ip — if'). The 
limit (9' if') — > (9 f), or equivalently — > 0, can be taken right away, it is not singular 
and it is just related to the spherical symmetry of the horizon, provided the coordinate y is 
cyclic. Then, since -Pz(l) = 1, 



1 1 t°° ■ 1 

G (x, x ) = — ■ : — — ■ ; — / dk x e X 
7r (1 + sm£j(l + sm£'J J-oc 

£(21 + 1) . = ■ e-i? ™«™~. (57) 

1=0 ^Hi)Wi(i') 

Let us start now with the point splitting procedure. We assume that the points x and 
x' are connected by a non null geodesic in such a way that they are at the same proper 
distance e from a third midpoint x. We parametrize the geodesic by its proper distance r 
and denote the end points by x and x' (it should not be confused with the third component 
of (£, 7], x, y)). Then we expand, 

e 2 

x» = ar"(r) = x^(r + e) =4 + e^ + -^ + ---, (58) 

a/" = x^(r') = x* 1 (f - e) = a# - ex? + - x£ , (59) 

where we have defined 



dx^ 

Xq = X>*(t), X^= — 

dr 



^ 

x 2 = 



d 2 x fl 



dr 2 



(60) 



i.e. the subindex of a given coordinate indicates the number of derivatives with respect to 
r at the point f. 

Using these expansions in e we can write /|, Wi(£")d£"* , [W / i(£)W / i(£')]~ 1//2 and the pref- 
actor [(1 + sin 0(1 + sin£')] -1 in Q57|) up to 0(e 4 ). Note also that x-af = 2 ex t + (e 3 /3) x 3 + 
(e 5 /60) a; 5 + 0(e 7 ). With these expressions we can expand the exponential term in fl57D in 
powers of e as 



-if*, Wi(€")d£"*+ik x (x~x') = e -ie(5i a;-<52 fc*) x Q _| \ 



where we preserve from the expansion the zero adiabatic terms because they will be useful 
as new integration variables. These terms have been included in the coefficients Si and S2 
as, 



e 4 „ _ e 2 e 4 



^ = 2£* + — £* + — <5 2 = 2xi + — x 3 + — a£. (61) 
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Since the number of derivatives d/d^* determine the adiabatic order, we introduce a new 
parameter T, which will indicate the adiabatic order, then at the end of calculations we will 
take T = 1. We will also denote, 



w = v l = v 1 (0, w = w 1 {0, w = 



dw^) 



Recall that from section 3 we know the function W\ up to adiabatic order four, thus from 
(1113) we have now, 



W = uj + T 



A, 



+ 



B- 7 



rp4 



A, 



+ 



Ba Ca 



Da 



OJ' 



uj 9 + uj 11 I 



where the coefficients u>i, A n , B n , C n , are given in (\L7[ ) and (p0|). 

We substitute these expansions in (|5T| ) and separate the different adiabatic terms by 
expanding in powers of the parameter T. The process is quite simple but tedious and the 
result is, 
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71 



l + sinf) 1 



2 (t*2 (3-sinQ cos 2 £ 



cos 2 £ 



x / ^ (2/ + 1) e" 

'- 00 1=0 



-ie(8i L0—S2 k x ) 



M 2 (l + sin£) 7 ^M(l + sin£) 5 / 
Aq + it A x + e 2 .A 2 + ie 3 A 3 + e 4 A 4 + 0(e 5 



+ 0(e 4 



(62) 



where the coefficients ^4 n are given in Appendix C. 

Let us now proceed to the integration in (0), evaluating first the summation over the 
discrete parameter I. It is convenient to transform the above summation into an integral by 



means of the Euler-Maclaurin summation formula [34|, 



/•« 1 1 1 

5>(0 = / TO-nn)]-- [F'(0)-F'(n)] + — [F"'(0) - F"'(n)} + - ■ • , 

(63) 

where F' = dF{l)/dl. In our case ra — > 00 and F^'(l) is given from (|62D as, F^(l) oc 
exp (— ie^icu). Note that in the limit n — > 00 with finite geodesic distance e, the phase in 
fW(n) is highly oscillating but the modulus is either aproaching to zero or is finite. In other 
words, the contribution at n — > 00 is negligible, this can be mathematically reproduced by 
introducing a small negative imaginary part in the geodesic splitting as, e — > e — z0 + , which 
allows us to regularize from the very beginning: lim^oo F^(n) = 0. Therefore we write, 

1 /-n 1 1 1 rn 

E *"(0 = y o ^ + 2 - u F ' (0) + 720 F "' (0) + " ' = Jo F(x) dx + ' 

(64) 

where the definition of the differential operator C should be clear from the last identity. 
Then we can substitute the summation in (|62|) for an integral over a formally continuous 
parameter /. This integral can be further simplified by using a new integration variable u, 
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uM~ l = lu = M -1 



kl + Vr 1/2 , ^ = a (0/(/ + l) + c(0, 



(65) 



where we define, see (|H|), a(£) = cos 2 £(l + sin£) 4 , c(£) = 2cos 2 £(l + sin£) 5 . Then 
2udu = a(£) (2/ + 1) eZZ, and, 



z(z + i; 



w 2 -f] 2 
*(0 



, ^ = A; 2 + c(e), 



where for simplicity we introduce a new dimensionless variable f2. With this one can directly 
substitute the derivatives d l Vi/d^* 1 = Vj which appear in the coefficients A n by, 



' M 2 
With all this (1621) takes the form, 



am 



(66) 



G^(x,x') = — 

7T 



1 



1 + sinO 2 

00 2-u 



* 2 ( 3 - sinf) cos 2 £ cos 2 £ 



M 2 (l + sin0 7 ?2 M(l + sinOV 1 



x 



dk x e ie&2h . 
oo Jn a(£) 



due -ie6 lU /M 



A + ie Ai + e 2 A + ie 3 A + e 4 A + 0(e 5 



+ 



(2Z + 1) (A + ie A + e 2 A 2 + ie 3 A 3 + e 4 A + 0(e 5 )) e"^ 1 ' 



1=0 



(67) 



where the functions A n are now written in powers of u 1 and Q. The integrals that we have 
in (|S7D are of the type,, 



21' 



30 /"OO O m 

dfci / e - ie{5lM/M - &fe} — du, 
oo w n 



(6? 



which can be easily integrated over u using the following identities, 



-ieS lU /M udu= _ 



M 2 

+ 1- 



-ie <5i Q/M 



e 2 8( e6 



-ie5 lU /M du 



iM 



-ie (5i Q/M 



1 U n 



.e5 l 



n-l 



H^jJ E i (-ie8 1 Q/M) + 

-ie Sl n/M 1 L_^V_ z £^y (n-2-fc)!; n > 0. 



(n-l)! fi r 



fc=0 



M 



where we have used a small negative imaginary part in e, i.e. e — i0 + , and where Ej(ix) 



is the integral-exponential function with imaginary argument [f4Q[ . The remaining integrals 
are, 



7--m _ 7- 
-'"0 — m 



e -ie( 5l n-* a fea )/M ^1 = ( e (m-l)) £ > o 
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T^=X m = e-^ n - 52k ^/ M Q m dk 1 = 0(e-( m+1 )) + J ?}z n( 2 x, 

J —oo I [1 T C/(£ JJ 



0(e) £ m > even 

m > odd, 



/oo 
e^^^E^-ie^fi/Mjfl" 1 * = 0( e - (m+1) ) + 0(e) £, m > even, 
-oo 

where i is a logarithmic term defined as i = 7 + ln(e/M), where 7 is Euler's constant. These 
integrals can be calculated using the indications of Appendix D. 

Many of the terms in ([)7|), however, are not relevant for our calculation because they give 
results of order beyond e 2 . In fact, note that the series Ao+ie A\+e 2 A2+ie 3 A^+e A Ai+0(e 5 ) 
has dimension of M" 1 , Q and u are dimensionless variables, T has dimensions of M _1 , and 
the geodesic coefficients d n x^/dr n have dimensions of M~^ a ~ l \ Now a simple dimensional 
analysis shows that the singularity a~ 2 in the Hadamard function is recovered with the single 
term Ao up to adiabatic order zero only, the logarithmic singularity, i.e. In a, is recovered 
also with the single term Ao, but now up to adiabatic order two, and the singularity a In a 
is recovered with the terms Ao, At, A2 and .A3 up to adiabatic order four. 

We can now substitute in (|67|) the values of the integrals given in Appendix D, the values 
for the functions a and c given in fl55p, the expression (pTI), for the coefficients 8\, 62, and 
the relation between the geodesic coefficients £* i x n in terms of £* and X\ (see Appendix E). 
Also we use the identities, see Appendix B, o = 2e 2 S, o^ 1 = 2ep' 1 , p^pp, = E, where a' 1 is the 
geodesic tangent vector in the midpoint x with modulus the proper distance on the geodesic 
between x and x' . In particular, a x = 2ex\. Finally, after a rather tedious calculation we 
obtain the expression for the Hadamard function, 

4 M 2 C 2 

G^(x, x') = +A + aB + Cxx a*a s + D s 



ira a 
where the normalization constant C, by comparison with ([45]), is necessarily C 2 = (16M 2 7r) _1 
and the explicit expressions for A, B, C ss and D gsss are given up to adiabatic order four 
in Appendix C. Note that, as it is expected on general grounds, there is no logarithmic 
singularity. In fact the coefficient for the logarithmic term is, 

2 c M 2 a 2 M 2 a 

3 a 3 a 3 3 a 2 ' 

which is identicaly zero as one can easily see by direct substitution of the value for the 
functions a and c, in flBSP- 

According to (|49|), the Hadamard function can be regularized using the elementary 
Hadamard solution PSf), which in our case reduces simply to, 



S{x,x')\ 



U'() 



4tt 2 (7 : 



since the coefficient of the logarithmic divergence, see fl47|) , exactly cancells up to order a 
for a Schwarzshild geometry (the interaction region is locally isometric to Schwarzshild). In 
fact, for our geometry and using 
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-oi - - a 2 a = — R aP \ R M a v - - — R a ^ 5 R a ^s <r = 0, 
1 2 180 M w 4 180 w 

Finally the regularized expression for G^(x,x') correctly normalized up to order e 2 is, 

X X X X 

G$(x, x')=A + aB + C^ a s a s + a ° ° - . (69) 

a 

From this expression we can directly read the regularized mean square field in the "in" 
vacuum state as (4> 2 ) = A. 

9 Expectation value of the stress-energy tensor 

To calculate the vacuum expectation value of the stress-energy tensor near the horizon we 
have to apply the differential operator ([0]) to (|69|). As we have already pointed out, this is 
not straightforward because we work with nonlocal quantities. Note first that the operator 
(H)) acts over bitensors which depend on the end points x and x', but the expression (|69|) 
for G$ depends on the midpoint x. It means, therefore, that we need to covariantly expand 
( jB9|) in terms of the endpoints x and x' . Consider the general expression, 

G$ (x, x') = A + a B + <x V + D &m , (70) 

where it is understood that A, B, C a p, are functions that depend on the endpoints x, 

x' but are evaluated at the midpoint x. Then, following the formalism of Appendix B, we 
can expand in the neighbourhood of 

Gg> (x, x') = A + a B + gl 4 C ^ ^ + < 4 % flJ D^s (71) 

where before expanding, a homogeneization of the indices by the parallel transport bivector 
g£, i.e, = cjp(T u , has been applied. The covariant expansions at x are given by, 

A = A - l -A, a a a + ~ A. a/3 a a a? + 0(e 3 ), B = B + 0(e), (72) 
2 o 

9a 4 = C aP + O(e), g% 4 g% g\ D^ s = D Q ^ S + O(e). (73) 

Now we can apply the differential operator to (|7T|). We will consider the two 
physically relevant cases of the minimal coupling (£ = 0) which should provide a good 
qualitative description for gravitons, and of the conformal coupling (£ = 1/6) which should 
provide a good qualitative description for photons (the use of £ as the coupling parameter 
should not be confused with the coordinate £ of the interaction region). If we introduce the 
operators, 

= V„/V„ + V„V„, C$ = V,V„ + V„V„, (74) 
the operator (H4T), for the minimal and conformal cases, is 
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-r,5=i/6 _ ±_ 
^ 12 



^ 4 



o a0 r(X) 



*/3 /'(I) 



a a afi £ (2) 



By the properties of the geodetic interval bivector a^ 1 given in the Appendix B, which can 
be also written as = 2ep M (p M £^ = S) we can prove the following identities, 



- lim, £g> a = -2 ^, lim, £« 



lim ; 4 2 )a^ = -4 5^^, 



C7 



- lim 6 2 } 



[IV 



a 



-% 5 { " & pV> + 64 6^ p% + 8 pVpV 



2 (a* 



a /3 7 (5 



where (• • •) is the usual symmetrization operator. Note that by a straightforward application 
of Synge's theorem (see Appendix B), the coincidence limits of Cfy and differ in a 
sign when they are applied over bitensors for which the first covariant derivative has null 
coincidence limit. But when they are applied over the biscalar A in (|72|) they coincide: 
lim x ^ x i Cfy A = ]hn x -* x > Cfy A = A ;fll/ /2. The application of the opperators ([74]) over quartic 
products of <7 M gives path dependent terms when the limit x —> x' is performed. Therefore 
the averaging discussed in Appendix B is required. Provided we work with a diagonal metric 
(as is the case), we only need the following averages, 



(pV> 



■sr 



(p x p x p x p 11 ) = - (g 1 



( P x p x p x p x p^) = ^ {g xx f 5^8 XV + ^ gT {g xx f ■ 

We now have all we need to evaluate the vacuum expectation value of the stress-energy 
tensor by application of the differential operator ( [4*4]) to the expression ( |7T| ) of G% (x, x'). In 
the orthonormal basis, 8q = M (1 + sin£) d£, 9\ = M(l + sin£) drj, # 2 — cos£(l + sin£) _1 dx, 
9 3 = (1 + sin£) cos rjdy, we obtain the following expectation values (TE,) for the minimal 
and conformal couplings 



x, x 



diag(pi, -pi, f>2, -Pi] 



(T B )C =l/6 = lim P J=l/6 G g) (X)X /) = diag(p; _ pj 3p) _ p) 

where p, p\ and p 2 are positive definite functions given by, 

cos" 4 £ 10771 



(75) 
(76) 



Pi 



P2 



256tt 2 M 4 2880 
cos~ 4 £ 3341 
' 256vr 2 M 4 9216 



1 + Ai cos 2 £ + 0(cos 4 «£) 
1 + A 2 cos 2 ^ + O(cos 4 
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p 



cos" 4 £ 6869 



(l + A cos 2 ^ + O(cos 4 



256tt 2 M 4 2880 

which are unbounded at the horizon (£ = 7r/2), and the approximate values for Ai, A2 and A 
are Ai = 1.474 • • •, A2 = 1.466 • • • and A = 1.469 • • •. From fl5lf) this is not the final result of 
the stress-energy tensor in region I, (T^j, since the trace anomaly term must be included. 
But this term depends on the spacetime curvature which is finite at the horizon, 



Q2 

64 tt 



9iw 



1 



- + O(cos 2 



diag (1, 



.240M 4 tt 2 64 

and it does not modify the dominant contributions of (|75j) or (ff^), which increase as cos _4 £ 
near the horizon. In the conformally coupled case the trace (Tjf)j is finite and can be 
obtained from the previous term. Inspection of (^5|)- (|76|) shows that for both couplings the 
weak energy condition is satisfied |43[], which means that the energy density is nonnegative 
for any observer, and that the strong energy condition is only satisfied for the conformal 
coupling. 



10 Conclusions 

We have evaluated the expectation value of the stress-energy tensor of a massless scalar 
field in the state, which corresponds to the physical vacuum state before the collision of 
the plane waves. This vacuum state has an analogue in the Schwarzchild black hole case 
as the empty state at large radius from the hole, the Boulware vacuum state, for which the 
vacuum expectation value of the stress-energy tensor also diverges at the horizon. In this 
case, however, it is argued that the Boulware vacuum state is not physical since it does not 
correspond to the vacuum of the gravitational collapse problem [f2B], |2"5| . 

Our results are the following. Before the collision of the plane waves (in region IV) 
(T Atw )/y = by the definition of our vacuum state, in the plane wave regions (regions II 
and III) we have found in section 7 that {T^ u )n = {T^ u )ni = 0, since the plane waves do 
not polarize the vacuum, and finally, in the interaction region (region I) (T^ v )i becomes 
unbounded at the regular Killing- Cauchy horizon (£ = vr/2) in the conformally and mini- 
mally coupled cases. In both cases the weak energy condition is satisfied, the rest energy 
density is positive and diverges as cos~ 4 £ and two of the principal pressures are negative 
and of the same order of magnitude of the energy density. The strong energy condition is 
satisfied for the conformal coupling, in this case {T£)i is finite but (T^) {T^J) 1 diverges at 
the horizon and we may use ref |41| on the stability of Cauchy horizons to argue that the 



horizon will aquire by backreaction a curvature singularity too. Thus contrary to the simple 
plane waves which do not polarize the vacuum [|1|, |T3|| , the nonlinear collision of these waves 
polarize the vacuum and the focusing effect that the waves exert produce at the focusing 
points an unbounded positive energy density. Therefore when the colliding waves produce 
a Killing-Cauchy horizon that horizon is unstable by vacuum polarization. 

In the more generic case when the wave collision produces a spacelike singularity it seems 
clear that the vacuum expectation value of the stress-energy tensor will also grow unbounded 
near the singularity. The reason is that such unboundness is, essentially, a consequence of 
the blueshift suffered by the mode solutions as they enter the plane wave regions, and it is 



easy to see |T6[ that any plane wave produces a similar blueshift on mode solutions. In view 
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of our results it seems very unlikely that the negative pressures associated to the quantum 
fields could prevent the formation of the singularity. 
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A Adiabatic expansion of "in" modes in region I 



With the parameters of section 5, defining p = Lg~ 1 (cos^) _1 and p\ = Lg~ 1 (cos rj)~ l and 
assuming that p < 1, pi < 1, we can expand ( [21~D up to adiabatic order 4 in powers of p and 
P\. The process requires to expand the WKB functions fll8|) and (|38|). These expansions 



are rather easy to perform and, in particular, the expansion for (38) allow us to fix the two 
arbitrary coefficients C\ and C 2 by requiring that ( |3"8"D corresponds asymptotically to the 
function ( p3|) which is, modulus a constant factor, an associated Legendre function. Using 
that the asymptotic expansion of an associated Legendre function P/™(sin?7), in the range 
of values, Lm(cosr]) _1 <C 1, is 



Pf m (sin77) = V 



I- 1 / 2 
y/n cos r] 



cos 



+ 0(1 



-3/2i 



we find that 



Ufa) = V 2 ^ + 1 _ 1 p »(;+m)7r/2 L-i J W^( V )d v * , C_ 1 y+m e i / ( V ) d v * 1 

With these results, we can divide each term (|3l|) or (j32"l) in four integrals such that the 
variation of the phase of their integrand functions are given by 



2L 2 k_ ~ y ' 2L 2 k_ 



A{U) ±A), 4^±w£ 4) M 



and the modulus of their integrand functions are given, respectively, by 

FM {^j- =F W*>(«)) + F 2 («), Fx(u) =F < + F 2 (u), 

F(u) (^g- =F W< 4) (r*)) + F 2 (u), Fl (u) =F W?>(«)) + F 2 («). 
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Note that when a sign ± is taken in a phase term, then in the modulus term we have the 
oposite sign =)=. Up to adiabatic order 4 all these functions, i.e. W± (u), Fi(u), F 2 (u) and 
A(u), can be expanded in powers of p on the boundaries between regions II and III with 
region I, where £ = ±77 and thus p = p±, and where we use for simplicity z = sinw. We get, 



W[ 4) (u) 



wYHu) 



Fo(u) 



A(u) 



L-{ 2 + /^(2-^) + g(l + ,r-^ 



+P q 
U 2 



<£_ 

64 



12 + \2z 2 + z 4 ) + q 2 



V 2 



16 
L- 1 [p 2 



(6 + 7z 2 \ 

( / V 2 



V 4 



2048 



256 
^ 4 



(l + z) 4 (30 - 24:z + 7z 2 ) + 

+o(L-y), 



U 2 



32 y J 2 



+ 



+p 4 



Q 



^ (1 + #(30 - 24, + 7z 2 ) - ^ (6 + 7z 2 )\ - 



V 1 



2048 



:i+z)* + 



o(L-y), 



= LM\2+p 2 



4 



(2-z 2 



32 



:i+z) 



U 2 



+ 



+p 

5U 2 
16 

Lp 3 M 



r -2 2 



l + ,) 4 (34- 24, + 5, 2 )- 



(2 + z 2 ))+^-(l + zf + 



5 f/ 4 V 2 U 4 



16 



128 



I + 



4096 

jz - ^ (1 + z) 3 (2 + z - z 2 ) + tf 2 z j + 0(Lp 5 ). 



0(Lp<), 
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:i + z) 4 + f/ 2 



Among the integrals in (^TJ)-(^), there is only one integral which contains stationary 
points and it is always in the term A^! in ( |31"D when either L _1 = L 2 k_ or L" 1 = L\k + , 
because the variation of the phase takes in this case its minimum value with respect to the 
modulus, and its contribution becomes dominant; consequently we can neglect the remaining 
integrals. The dominant integrals connecting the "in" modes between regions II and I can 
be written as, 



du Fi(u) 



+W {4 \u)) e^ ' du ( A(M)/(2i2M_W/ - 4){M) ) 



2Lok- 



which can be written in powers of p as (0). The dominant integrals connecting the "in" 
modes between regions III and I can be written in the same way but changing L 2 k_ for 
Lik+. 
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B Bitensor algebra 



A bitensor with n + m indices, a x , a 2 , ■ ■ ■ , a n , (3 1 , (3 2 , ■ ■ ■ , (3 m , T aua2 ^.. tanl3ul32 ^.. tl3m , is an ob- 
ject which transforms, under coordinate changes, as the product of two ordinary tensors, 
A ai a 2 -a n an d Bpp p m , each at a different spacetime point |35|, f|2|, |43|. The standard ten- 
sorial operations can be straightforwardly extended to bitensors with the simple precaution 
of contracting only indices refering to the same spacetime point, and noting that when a 
covariant derivative is performed over one of the points, the indices referring to the other 
point are irrelevant. 

If we have an arbitrary bitensor evaluated at two close spacetime points x and x' it is 
possible to expand it in a covariant way at one of the two points. To achieve this we make use 
of the geodetic interval, s(x, x'), a basic object in geodesic theory, which describes the proper 
distance between the spacetime points x and x' along a non null geodesic connecting them. 
Let us restrict ourselves to a neighbourhood of x (or x') where the geodesies emanating from 
x do not cross (mathematicaly this means that the unique geodesic joining x and x' is an 
extremal curve jHJ), then s(x,x') is single- valued. The basic properties of s(x,x') are 

g^s^s,* = g^'sysy = ±1, (77) 

together with the symmetry s(x,x') = s(x',x) and \im x ^ x / s(x,x') = 0. 

These ensure that the bitensor g fiV (x)s tfM (x, x') is a tangent vector to the geodesic at the 
point x and a scalar at x' . Note that g tiu (x)s^(x,x') is an orthogonal vector to the surface 
s(x, x') = constant, and this surface is orthogonal to the geodesies which emanate from x' and 
therefore g tJ,v (x)s t ^(x,x > ) is a tangent vector to the geodesic that crosses s(x, x') = constant 
at the point x. Of course, if x is fixed and we draw the geometrical locus of the points x' at 
the same proper distance s(x',x), from x an equivalent description is found. If we take the 
signature (+, — , — , — ) then the positive sign in ( |77| ) refers to timelike separated points, and 
the negative to spacelike separated points. When s(x,x') = 0, then x and x' lie on a light 
cone. 

Instead of s(x,x') it is convenient to define the geodetic biscalar, o~(x,x'), which unlike 
s(x,x') has no branchpoints, as 

a = S - S, (78) 

where we introduce the notation £ = ±. For the signature (+, — , — , — ) the plus sign is for 
timelike separated points and the negative sign for spacelike separated points. Let us define 
o~ [i = cr M , then the following relations are trivially satisfied, 

~^<7^ = ^ V V<V = *» ( 79 ) 

together with the coincidence limits \im x ^ x r o = 0, and lim x ^ x i a tlM = lim x ^ x r o~^> = 0. 

Given an arbitrary bitensor A^x, x'), with both indices referring to the same point x, 
it is possible to covariantly expand this tensor in the neighbourhood of x, as 

A a(3 (x, x') = A aP {x) + A 7 Q/3 (x) cx 7 + ~ A* 6 a p(x)a 7 a s H , (80) 
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where the coefficients A a p{x), A t a p{x) : • • • are ordinary local tensors at the point x. These 
coefficients are easily calculated, provided that the successive covariant derivatives of the 
geodetic biscalar a are known. For this we derive ( |T9"P recursively: 

0> = g a ^<7 a a-/3^, cr.^ u = g a/3 (a- au a-i3^ + a a cr ]l3fll ,) , (81) 
Then it is straightforward to find, 

where the notation introduced by Synge to indicate the coincidence limit x — > x' is used, 
i.e, [T aijai> ... itW 3 lAj ... i/3 J = \xm x ^ x >T aiiO , li ... tOln li j 3u ... i p m (x,x'). In particular, it is possible to 
covariantly expand the successive covariant derivatives of the geodetic biscalar a as, 



) + 0(s). 

Sometimes, for the sake of symmetry, it is convenient to covariantly expand the bitensor 
in terms of the midpoint x of the geodesic connecting x and x' . For simplicity let us take a 
biscalar G(x,x'). Then a midpoint covariant expansion reads, 

G(x, x) = A(x) + A„(x) a" + A {fJtu) (x) a" a v + A (liur) (x) ^ a u a T + ■ ■ ■ , 

where (• • ■) index symmetrization. The relation between the expansion coefficients A(x) and 
the biscalar G(x, x') is easily found to be, 

A = [G] , A, = [G.J - \ [G^ , 2! A, u = [G.^] - [G.^ u + \ [G]^ , 

3 3 1 

3! A^ UT = [G-^ VT ] - - [G. )IXU ]. T + - [G.J[. VT - - [G] ;iJa/T , 

where the notation = means that the equality is modulus terms which vanish under sym- 
metrizations. 

The covariant expansion of a bitensor with indices referring to different points is slightly 
more sophisticated because it is necessary to homogenize the indices before the expansion. 
To do this we use the parallel transport bivector, g^,(x,x'), which is defined by analogy to 
the parallel transport of ordinary vectors on a manifold. Given an ordinary vector A^{x) at 
the point x it is possible to parallel transport it at the point x' on a non-null geodesic, say 
A 1 - 1 (x'), by a linear operator g^(x, x') as A 11 (x') = gfi (x, x') A^(x), provided that g%{x, x') 
is covariantly conserved along the geodesic, i.e. cr v '(x, x') g^(x, x').y = 0, and it satisfies the 
trivial boundary condition, g£ = 5£ . In exactly the same way, we can define a parallel 
transport for bitensors as g^, (x, x') if it verifies the previous differential condition not only 
at the point x' but also at the point x, i.e, 
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a"(x,x') g^(x,x'). v = 0, 



and the same boundary condition 



The parallel transport bivector, defined in this 



a way, is unique because if we fix the point x one can integrate the second equation in (|83|) 
along each geodesic which emanates from x up to x' with the initial condition 



6" 



Reciprocally, if we fix x', the first equation in (|83| ) can be integrated from x' down to x with 



the same initial condition. This reciprocity can be expressed as, g p ,(x, x') 
the geometric interpretation of g^,(x,x') it is straightforward to see that: 



9u\ x , x ) 



From 



9u 9 V > S/jV — g^vi 9a' 9 V < 9nv — 9^'u 1 , 



9^i <V 



-a 



9u 



-<v> 



^ 9u' 



Thus given a non-homogeneous bitensor A^ v >(x, x'), we first homogenize it as, x') = 

(x,x , )A liV i(x,x r ), and then we covariantly expand A^x, x') in the same way as flSOD . 
An interesting result in the algebra of bitensors, usually known as Synge's theorem 
45j, is that given an arbitrary bitensor T aia2 ... an prp> p^, whose coincidence limit and the 



coincidence limit of its first covariant derivatives exist, then the following identity is satisfied 



42]: 



- T aia2 ... an p' i p' 2 ...p lm . fJi + T aia2 ... oln p> i p> 2 ...p, 



54) 



This is very useful when the calculations involve a lot of primed indices and we know the 
coincidence limits of the nonprimed quantities. A trivial application of this result implies 
that [<j^\ = 9 ufi 1 • 

When one is performing the covariant expansion of Green functions in preparation for 
the calculation of the regularized expectation value of the stress-energy tensor, ambiguities 
in the coincidence limits appear when one of the points, say x, is held fixed and an arbitrary 
point x' is allowed to approach x. These ambiguities are due to the different paths that x' 
may follow, therefore some type of averaging is required. The most elementary averaging 
is called four dimensional hyperspherical averaging and it consists in giving the same 
weight to all the geodesies which emanate from x as follows. First one analytically continues 
to an Euclidean metric the components of the tangent vectors to the geodesies which emanate 
from x. Second, one averages over a 4-sphere, and third the results are continued back to 
the original metric. It is not very complicated to find the following averaging formulae which 
are useful in this paper, 



1 1 5 

(c^ov) = ^V) (o- v o- T a p ) = -a 2 g^ u g Tp) , (a pO v o T o p o ^ = - a 3 g^g^g^)) 



Note that from the symmetry of the averaging procedure, the average of an odd number of 
components cr M vanishes identically. 
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C Useful adiabatic expansions 

The adiabatic coefficients A n in formula fl62|) are given by, 



An 



A, 



1 

- + 



5 V t 



1 V, 



32 uj 7 8 



/n^V^_231lf^ 21 
I 2048 0^ ~ 256 cj 11 + 128 



16 cu 6 4 w 5 / ^ 



7 W 



32 a; 9 



/u^J^_2jttlf^ 21 V£ 7 ^ 
I 1024 ^2 ~ 128 u 11 + 64 U + 16 





A. 



IM r + T2 (iv? 



35 V, 3 




4 cj 3 

1 



6 



*2 



128 cj 9 



525 V^ 4 



ii 



1^4 _ - | , 
16 uj 1 16 1 Vs - 



512 u 



+ 



217 V^ 2 VJ 



11 V? 



13 v,v, 



128 uj 9 32 cj 7 32 cj 7 



+ 



J_Vi 
16 oj 5 



*2 




1 M 

2^J 4l42 
35 

64 a; 8 



I* 

5 



96 w 6 24 w 2 J ?3 I 12 cu 4 6 cj 2 



+ 



1 V, 



ijj' 



SI S2 



7 Wz 



2048 ^2 _ 256 lu w + 128 ^ 96 ^7 



96 cj 6 



35 V? 63 IfV, IV, 2 llWi.l^,, 



*3 



64 uj 10 64 



4 u e 48 u; 6 



24 u; 4 



+T 3 

+T 4 
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1536 to- 



315 V 4 245 W, 5 V? 



1 Yl _ 1 M r 2 + ( 1H _ 1 H W 

32 cj 5 16 w 3 / ?2 1 12 cu 5 12 w 3 J ? 



15 Vf 9 VM 
"32 cj 7 16 uj 5 



1024 w 




+ 



15 ViVi 



1 V/ 



-*2 



512 cu 9 64 cj 7 128 u 



175 VJ 4 217 Vfg _ 11 V? 
512 ^ + 384 uj 9 ~ 96 To 



, .11 ooa . .Q , .7 gg ^7 + A8 ' ^ ^ 3 
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35 V* 37 Vf Vz 13 V? 1 VjVi 1 V, , 4 
192 cj 9 96 96 cj 5 12 cu 5 48 cj 3 1 SJ 



T 



12 



2 ?3 



24 cu 2 



The coefficients for the Hadamard function in using 2; = sin^ are: 



16MVA 
16M 4 7r 2 i? 



(19136 + 1523^ + 99^ 2 - 950^ 3 - 422^ 4 + 3z 5 + 3z t 



:i+*r 

14400 



;-1093271 + 1317560z - 275260z 2 - 187012z 3 + 37110z 4 - 



4908z 5 - 23962 b + 772^ 7 + 57z 8 - 12z 9 )(l + z 
2{4z - 3)(1 + z)-\l - In [(1 + z) 2~ 1/3 



-6 ' 



115200 



+ 



16MVG 



-1474763 + 1403169z - 366148z 2 - 39568z 3 + I86860 4 - 774^ 5 

z 2 



196z 6 + 56^ 7 + 21z 8 - 3z 9 )(l + z 
2(1 + z)' 6 In Ul + z) 2~ 1/3 



57600 



+ 



16M 4 7r 2 £> 



'-3-Az + 2z 2 y(l + z) 



\-8 



'1-Z 



2\-l 
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D Some useful integrals. 



Let us define 



/oo 
dk tt n e~ ie( - Sl Q - h k \ tt = Vk^TX; A > 0. 
-00 



15) 



For the particular value n = — 1 this integral can be easily solved. Performing the following 
change of variable, 



t = t(k) =5 1 n-5 2 k = 5 1 Vk 2 + X - 5 2 k, (86) 

we must consider separately the two possibilities \5i\ > \8 2 \ or \5i\ < \8 2 \. For \5i\ > \8 2 \ and 
without loss of generality we may take 8%, 5 2 > 0, and invert the change ( |36"D as, 



k = k(t) = 7 



8 2 t±8 1 Jt 2 -(3 



J7) 



where we use two new variables 7 and j3, defined as, 



7- 1 = 51 - 51 (3 2 = Xj-\ (88) 

Note that k(t), in (p7|) , is a bivalued function and we have to be careful in changing the 
integration limits. Since \5i\ > \8 2 \, then \im.k_ > ± 00 t(k) = +00, and this means that (BH) has 
an absolute minimum at the point (+7/3^2, /?)■ With this it is easy to see that we have 
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to take as inverse function of flSlf ) the function ( |5TD with the plus sign to the right of the 
minimum t = (3 and with the minus sign to the left of the minimum. Therefore we can split 
the integral (BH), for n = — 1, in two parts at each side of t = j3, i.e, 



J_ = P^ 2 * e -ie(6i n-6 2 k) + f°° ^ e -«(5i fl-S 2 k) 
J-00 f2 J^BSo. Q 



(89) 



f" 00 dk 

The change of variables (|86|) gives Vl~ 1 dk = ±(t 2 — P 2 )~ 1 ^ 2 dt, with the minus sign for the 
first integral in (|89|) and the plus sign for the second integral. Then the integration can be 



easily performed in terms of a zero order Bessel function 40 , as 



-ie t 



= dt = 2K (ie/3). 



(90) 



//3 v 7 ^? 2 

For the case \8\\ < \5 2 \, we define the parameter 7, as, 7 _1 = 5\ — Sf, and the inverse 



function of (B6|) is 



= fc(t) = -7 



(91) 



which is again a bivalued function. Now, however, the function ([86]) has no extrema, and it 
is easy to see that we can take a single inverse everywhere as, 



k(t) = -7 



S 2 t-5 1 Jt 2 + f3 2 



dk 



dt 



so that 

as, 



for n — — 1, can also be integrated in terms of a zero order Bessel function |40 



00 dfc 



-ie(Si fl— S2 k) 



-ie t 



dt 



cos (et) 



n ~ J-00 jWTW h ^t 2 + p 2 

Finally putting together the results (|90|) and (0), we can write 

J_ 1 = 2K (ie(3), 



dt2 K (e (3) 



with the parameter (3 given by, 



(3 2 = A7- 1 = A (51 - 51) . 



(92) 



(93) 



(94) 



The integrals X n with n > are divergent unless we adopt the standard prescription of 
taking e — ► e — z0 + for |5jJ > |5 2 | and e — > e + i0 + for |d>i| < | <5a | - In that case all these 
integrals can be easily related to Z_i by means of the following recursion formula, 



T 



i dl n - 



then Z„ 



n+l ffn+1 



(95) 



Since Z_i is essentially the Bessel function K and the derivatives of the Bessel functions can 
be written also in terms of Bessel functions ||40|| , then all Z n with n > can be expressed as 
a linear combination of Bessel functions. For example, using ([53) for /3 and recalling that 
K' (z) = dK (z)/dz = -Kx(z), 
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The integrals X n with n < —2 are finite in the limit e — ► and they can also be recursively 
calculated by means of Z_j and the following recursion relation for n < — 1: 

-QX=2 X ^- l — ln - 1 - (96) 
This allows us to calculate all the Z n with n < — 1 only if we know at least another integral 
besides Z_i. Fortunately it is not very difficult to evaluate Z_2 from the recursion relation 
(|9~5|) and by a simple integration, 



Z_ 2 (5i) =-ie I 1 dS[ J_i(<5j) +X- 2 (0), 

JO 



where the integration constant Z_2(0), which corresponds to the value of Z_2 for the special 
value 5i = 0, can be calculated directly from (|85|) ffl], then 



X_ 2 (5 1 ) = -i2e J*' ^K (ze/5(^)) + ^el eAl/2& l. (97) 
Let us denote, 

/oo 
dkn 2n Ei(-ie5 1 tt)e ie52k , Q = Vk^TX; A > 0, (98) 
-oo 

where Ei(zx), with x real, is the integral-exponential function with complex argument f40| . 
We can easily calculate X£ since by the properties of Ei(zx), it is not difficult to relate X£ 
with Z_2, 

dI£ ° - 1 X (W) 
-d\-2 T - 2 - (99) 

This can be easily integrated as 

l£o (A) = ]- f X d\'X^{X')+X£ (0), (100) 
2 Jo 

where, as in the case for Z_ 2 , the arbitrary integration constant Z£ (0) can be directly 
calculated from ( Pq) |40[ , for the particular value A = 0. The final result is, 



In 



(S 1 + 6 2 Yj , for S 2 ^Q. 



IS (X)= 1 - f X d\'X_ 2 (X') + \ e5 2 -"^-^' 'J' ' (1Q1) 



2 Jo ' ~, for 5 2 = 0. 

I ed 1 

The integrals Z£ 2 n with n > diverge unless we adopt the standard prescription e — > 
e — i0 + for |<5i| > |5 2 | and e — >■ e + i0 + for < |5 2 | then they can be calculated with the 
known values of X n and X£ Q by means of the following expression, 



<9<5? J-OO \ 
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which can be integrated by parts, using the properties of the Ei(ix), to obtain the following 
recursion relation 

IS n+ 2 = J" [In + (n + 2) ISo] . (102) 

E Geodesic coefficients £* and X{ 

We start with either a timelike or spacelike geodesic connecting the points x', x and x, which 
can be written in terms of the proper geodesic distance r as, 

x » = x m ( f + e ) j ^ = ^ ( f ) ; = ^ (f - e) , 

and define at the midpoint x the parameters £* = d n ^* /dr n \ f and x n = d n x/dr n \ f , which 
determine the geodesic. From the metric in the interaction region keeping 77 and y fixed, by 
the simmitry of the problem, we get 

^(ttS!)^* 2 -^)- (io3) 

If we parametrize the geodesic with its proper distance, then (|103|) gives the following relation 
between the coefficients £* and Xi, 

t? -£ = (104) 

1 — sm 4 

with the plus sign for timelike geodesies and the minus sign for spacelike geodesies. Since 
the coordinate x does not appear in (|103|) , it is easy to see that, 

dx_ ( l + sin A (d£\ 2 _ ( I + sin A 2 1 + sine 

rfr" Ml-sind' [dr "Ml-sinJ l-sinf 1 j 



Equations ( |105|) and the relation between £ and £*, d£* = d£M (1 + sin£) 2 / cos£, allow us 



to obtain the derivatives d n ^*/dr n and d n x/dr n , which when evaluated at r = r produce 
the geodesic coefficients we need. 
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Fig. 1: The four spacetime regions and their boundaries are represented in the coordinate 
system (u, v), which does not preserve the causal structure. Region IV is the flat region before 
the arrival of the waves, regions II and III represent the plane waves and region I is the interaction 
region. The boundaries C and £ ', which are spacetime lines, must be identified with the folding 
singularities V and V of the horizon u + v = n/2 (there is no spacetime beyond these lines). 



